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SECTION 1 


PURPOSE 


5° The purpose of this report is to describe an electromagnetic propagation conditions 
climatological database. a subset of Which is employed by the Historical Electromagnetic 
Propagation Conditions Summary Function (HEPC) within the Tactical Environmental 
. Support System (TESS) (reference 1). 


L& 


BACKGROUND 


The HEPC function generates a climatological description of the refractive condi- 
tions for a user-specified geographic location. The climatological description is composed 
of five parts: - 

5 eA < * . : 
® a percent occurrence of enhanced surface-to-surfave radar detection, electronic 
support measures (EMS) intercept and communications range 
f 
a surface-based duct summary , 
an elevated duct summary , 


an evaporation duct histogram > +. 


a genera] meteorology summary. 


Figure | illustrates the description. 


The statistics displayed by the HEPC function are derived from two meteorological 
f databases, the Radiosonde Data Analysis I] assembled by the GTE Sylvania Corporation and 
the Duct63 assembled by the National Climatic Data Center. Prior to the discussion of the 
HEPC database construction, it is appropriate to describe these two sources of data. 
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HISTORICAL PROPAGATION CONDITIONS SUMMARY 


Specified location: 20.00 N 430.00 W MS = 86 (° = Insufficient data) 
Radiosonde source. 30.00 N 140.00 W MS = 123 WMO aumner = 4YN 

Coastal station: F1XEO SHIP, NORTH PACIFIC OCEAN Hgt = 12m 
Surface observation source: 25.00N 125.00 W MS= 85 


Percent occurrence of enhanced surtace-to-surface radar/esm/com ranges: 


Average Feb 
nite 


Surface based duct summary 


Average Feb 
Parameter day nite d&n day nite 


Percent occurrence 28 3 


Avg thickness Km 
Avg trap freq GH? 
Avg tyr grad -N/Km 


Mar 


d&n 


d&n day nie d&n 


16 
09 
96 

220 


Elevated duct summary : 


Average 
Paramete: day nite d&n 


Percent occurrence 
Avg duct top Km 


Avg thickness Km 


Avg trap freq GHz 
Avg lyr grad -N/Km 
Avg lyr base Km 


Evaporation duct histogram in percent occurrence. 


Average 
Percent Occurrence nite 


Oto 4 Meters 
410 8 Meters 
8B to 12 Meters 


12 to 16 Meters 
16 to 20 Meters 
20 to 24 Meters 


24 to 28 Meters 
28 10 32 Meters 
32 to 36 Meters 


36 to 40 Meters 
above 40 Meters 
Mean height Meters 
= of observations 


General meteoro'cay summary” 


Average 
Parameter day nite d&n day nite 


* Accepted soundings 120 120 115 112 
©. occur EL&SB acts 

% occur 2+ El dets 

Avg station N 

Avg station -N/Kin 

Avg sfc wind m/s 


Figure 1. HEPC Function output for metric units. 
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SECTION 2 


GTE SYLVANIA RADIOSONDE 
DATA ANALYSIS I] DESCRIPTION 


General 


GTE Sylvania, under contract by the Department of Defense (DoD), conducted a 
large scale analysis of approximately three million worldwide radiosonde soundings from 
1966 to 1969 and 1973 to 1974. Numerous statistics of tropospheric ducts and super- 
refractive layers (SRLRs) were compiled. The basic pressure, temperature, and dew-point 
temperature profiles for each of the soundings Were expanded and converted into a refrac- 
tive index profile with 2 millibar increments. A search of each expanded profile was con- 
ducted to locate occurtences of gradients less than -157 N;Km Cin the case of ducts) or less 
than -100 N/Kin (in the case of SRLRs). Monthly estimates of the probability of occur- 
rence of ducts and SRLRs were then computed for each of the 92] radiosonde stations. 
Characteristics such as height. thickness, intensity, and minimum trapping frequency were 
also computed. 

The radiosonde Long-A tape. a truncated version of the Long Plus 2 tape data set 
represents the final output for the GTE Sylvania project (reference 2). The Long-A tape 
consists of 1842 records containing tropospheric ducting characteristics and weather-related 
parameters for 921 radiosonde stations distributed worldwide. Each set of two consecutive 
records supplies all the data for a particular station. The records are formatted in ASCII 
code with a length of 7104 bytes each. Appendix A lists the format for records ] and 2 for 
an Observation station. Appendix B is a listing of the included radiosonde stations. 


Term Definitions 

Layers are defined as surface (S), elevated-surface (ES), or elevated (E). Figures 2, 
3, and 4 illustrate the M-unit versus height relationships for these layers. 

The following is a description of each item contained within the two records. 
Appendix C contains a listing of the data as it appears on the GYE/Sylvania Long A data 
tape for station 4YN, fixed ship, North Pacific. 


Station Number. World Meteorological Organization (WMO) block and station 
number. 


Station Name. Station name by city and country together with a numerical indica- 
tor of coastal (1) or inland (0). 


Record Sequence Number. Either record | or 2. 
Latitude; Longitude. Station latitude and longitude to hundredth of a degree. 
Positive values represent north (west) and negative values represent south 


(east) latitude (longitudes). 


Elevation. Station elevation above mean sea level in meters. 


Local Time. Local time corresponding to midnight Greenwich Mean Time. 
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Figure 2. M-unit versus height profile for a surface duct. 
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Figure 3. M-unit versus height prufile for an elevated surface duct. 
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Figure 4. M-unit versus height profile for an elevated duct. 
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Local Time. Local time corresponding to noon Greenwich Mean Time. 


Station Elevation Reliability. Errors in station elevation contribute directly to 
etrors in duct heights and slopes. Various kinds of discrepancies were noted in 
the elevation of each radiosonde station. For example, during the measure- 
ment period, the station may have changed locations and thus elevations. 
These discrepancies have been compensated by a correction scheme. Values 
of Ltmost reliable) to 7 (least reliable) have been applied to the results of this 
compensation, 


Instrument Type. A numerical designation of the two most popular types of instmu- 
ment used during the list 2 yeurs together with the percentage of the data 
dunng those 2 years for which those instruments were used. Table 1 provides 
the definition: of the instrument designators, 


Table). Instrument types definition. 


Indicator Detinition 
6) Uninown 
) USWB exposed thermister -403MC21680MC R/S 
2 Australian-Variable AM R/S 
3 Bendix AN/AM1-4/46/4D; 12; GMD-1 TMQ 3/TYPE 430 
4 French -- Measural 
$ Canadian - Sangamo 
6 Indian Chron — AM R/S 1680 MHz 
a Indian FAN -- R/S 
Es) Russian A 22 Malahut/RKZ-2 
2) Chinese 
10 West German — GRAW M60 
11 Pakistani FM RS 403 MHz 
12 Finnish -- Vaisala 
)3 Japanese - Code sending 
14 Malay s:a Asior 403 
15 Hahan - Autovox IA‘ AMT 
16 East German — Freiberg RKS2 
1? British - KEW MK IB 
18 Dim Ty pe R.V. 4 
19 Omera Decea 
20 PEU Swiss Institute 
aa MARS 1K 


Accepted Soundings. The number of soundings tor which an expanded retractive 
index profile was computed and a search for ducts and SRLR was conducted. 


Accepted Soundings With Surface. Number and percent of soundings for which 
complete surface data were available. 


Accepted Soundings With Signal Level. Ducts are separated into those whose 
minimum trapping frequency exceeds 10,000 MHz and those with frequency 
less than 10.000 MHz. Since measurement errors and approximations inherent 
in the computational procedures are probably responsible for the occurrence 
of the thin, high frequency ducts, only the ducts with minimum trapping 
frequency less than 10.000 MHz are included in the statistics. 


Accepted Soundings With 1000-mb Winds. Number and percent of soundings which 
contain both wind speed and direction for the 1000-mb level. 


Accepted Soundings With 850-mb Winds. Number of soundings which contain both 
wind speed and direction for the 850 mb level. 


Average Mandatory Levels. Average number of mandatcry levels up to the 100 mb 
level contained within the soundings from the last 2-years data. 


Averaye Significant Levels. Average number of significant levels between the surface 
(but not including the surface) and the S00-mb level contained within the 


soundings from the last 2-years data. 


Pressure. Median value of the attnespheric pressure in millibars at the surface of the 
reporting station. 


Temperature. Median value of the atnospheric temperature in degrces Celsius at 
the surtace of the reporting station. 


Dew Point. Median value of the atmospheric dew point temperature in degrees 
Celsius at the surface of the reporting station. 


NS. Median value of the refractive index at the surface of the reporting station, 


NW. Median value of the “wet component of the refractive index at the surface of 
the reporting station. 


dN/dZ (xx). Median value of the refractive index gradient over the first xx meters 
above the surface at the reporting station. 


Wind Direction (xx). Estimate of the wind direction distributed over a circle at the 
xX pressure level. 


Wind Speed (xx). Median wind speed in knots at the xx pressure level. 


Number of Ducts (5 Years). Number of ducts occurring overt a 5-year period of the 
data, There is a separate entry for E, ES. and S ducts. 


Number of Soundings With Ducts. Number of soundings containing at least one 
duct of the type E and a combination of S and ES. 


Number of Soundings With SRLRs. Number of soundings containing at least one 
refractive layer of the type E and a combination of S and ES. 
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Percent Soundings With Ducts. Percent of soundings with ducts of type C, Es, or S. 
Statistics are provided for various time intervals and duct occurrences. 


Probability >= 2 Elevated Ducts. Probability of 2 or more EF ducts occurring 
sinultaneously. 


Probability S and E Ducts. Probability of S and E duct occurring simultaneously. 


Median Maximum Height. Median maximum height in meters of all ducting 
Occurrences. 


Coupling Factor. The difterence between the M-unit value at the surfac2 and the 
top of the duct 


Bottom Height. The height in meters from the surface to the bottom of the duct. 


Optimun Coupling Height. The height in meters from the surface at which the 
M-unit gradient becomes just trapping (inflection point). 


Top Height. The height in meters from the surface to the top of the duct. ForS 
duets, this is equivalent to the optimum coupling height. 


Thickness. The distance in meters from the bottom to the top of the duct. 


Intensity. The ditference berween the M-unit value at the inflection point in the 
profile and the top of the duct. 


Frequency. Median, standard deviation, and mode of the maximum electro- 
magneto frequency trapped within the trapping layer. 


Average Gradient. Median of the M-unit gradient within the duct. 


Delta M. Phe difference between the M-unit value at the surface (for S ducts) or at 
the inflection point (for ES and E ducts) and the top of the duct. 


Upper Slope. The M-unit gradient over the distance from the inflection point to 
the top of the duet, 


Lower Slope. The M-eunit gradient over the distance from the bottom of the duct to 
the anflection point, 


Surface to Inflection Point. The heiht in meters from the surface to the inflection 
pomtan thy refractive index profile. 


Surtace to Bottom, The height in meters trom the surface to the bottem of the duct. 
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deviation is the unbiased sample standard deviation of the value of the duct parameter. The 
CI (mean) gives the width of the 50-percent confidence interval for the mean, assuming that 
the sample was taken from a normal (Gaussian) distribution. The confidence interval width 
for the median will, in general, be wider than the confidence interval for the mean, since the 
median confidence interval is based only on the relative weak assumption of symmetry of 

the underlying distribution. For the case of wind direction, the 25-and 75-percent confidence 
interval limits are based upon symmetry about the estimated wind direction for the xx 
pressure level. 


Statistics 

The mean is the average of the values of the given duct pararneter. The standard 
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Additional Considerations 
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a. Where monthly variables are found in groups of 13, the order of the data is for 
January through December and annual total. 


ete 


b. Where monthly numbers and percentages appear in groups of 7818 values, 
the order is 


ce. Where groups of 3818 values are given for medians and confidence intervals, 
the order is 


bt 
ne I. number then percentage, alternating for months January through 
ry December and annual total. 

ae : es | 
a 2. same as (1) above but for time period 0000 GMT only 
aa 2 2 z { 
~ 3. sume as (1) above but for time period 1200 GMT only | 
wet | 
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1. the lower confidence intervals for January through December and 
annual totals. 


2. the medians for 12 months and annual totals 


a) 


the upper confidence intervals for 12 months and annual totals. 
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SECTION 3 


Ss: 


DUCT63 DATABASE DESCRIPTION 


sh 


fe General 
al The National Climatic Data Center, Asheville. North Carolina, under contract from 
mh : the Naval Ocean Systems Center (NOSC), produced a subset analysis of its Standard Tape 
| ‘ Deck Family bt (STD-11) database. The STD-11 database consists of over 150 years of 
Tp worldwide surface meteorological observations. These observations were assembled fram 
we » ship logs, ship weather reporting forms, published ship observations, automatic buoys, 
aM teletvpe reports, and card decks purchased from foreign meteorological services. 
an The subset analysis, known as the DUCT63, covers 293 Marsden squares and spans 
es 15 vears of surface observations. A Marsden square is a region of the earth’s surface defined 
i by avrid of 10 degrees latitude by 10 degrees longitude and is assigned a unique identifica- 
a tion nuraber. Figure S shows the location and the numerical assignment of all Marsden 
wn squares. For example, Marsden square o is defined as the region bounded by the prime 
a meridian to 10 degrees west longitude and from the equator to 10 degrees north latitude. 
a Not al! of the 648 possible Marsden squares are included in the DUCT63 analysis-for two 
r reasons. First, the analysis is specifically concerned with the maritime environment. : 
“Te Marsden squares not containing a region of ocean are excluded from the data. Second, a 
vy requirement of atleast 100 valid observations per month was imposed to reduce the effects 
ie of any spurious meteorological measurements on the distributions. 
£0 Figure 35 shows the location of the Marsden squares contained within the DUCT63 
oS analysis as the region enclosed by the heavy border. 
‘ae DUCT63 Database Elements 
n The DUCT63 analysis contains distributions of meteorological quantities and 
p surface-to-surface attenuation rates for frequencies of 35 and 94 GHz. These distributions 
*) ure expressed as either a probability or a percentage of time that the quantity is observed 
he within a specified range. The distributions include diurnal effects where day categories 
a imply a positive solar angle within the Marsden square at the time of the meteorological 
A Observation, Night categories are times of observation between one hour after the local 
; oud sunset und one heur before the local sunrise. Observations taken in the interval between 
nut the day and night categories are excluded from the data set. 
: There are two forms of distributions: a probability (or percent of time) distribution 
MY): fora specific quantity (i.¢.. wind speed) and a joint probability distribution of two quantities. 
ra The latter form, also called a cross distribution, is specifically desizned for use by NOSC in 
x4 : its research ciforts. 
x The following quantities are distributed within the DUCTO6%5 analysis and appear 
6. within the data in the order listed. Appendix D contains a complete fisting of the data for 
‘ " Marsden square number 8§, 
, 
: ’, l Paulas evaporation duct height (meters). 
Ms 2. Paulus evaporation duct height crossed with Jeske duct height. 


3. Wind speed (meters por second). 
4 Absoiute humidity (siums per cubic meter). 
et 
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Figure 5. Marsden square numbering system for the worid. 


>. Modified air-sea temperature difference (deg C). 

6. Rain rate (scaled millimeters per hour). 

7. Attenuation rate (Gaseous) at 35 GHz (scaled dB/km). 
8. Attenuation rate (Gaseous) at 94 GHz (scaled dB/km). 
9. Attenuation rate (total) at 35 GHz (scaled dB/km). 
10. Attenuation rate (total) at 94 GHz (scaled dB/km). 


11. Paulus evaporation duct height crossed with wind speed. 


Each Marsden square’s data are contained within a record of 55242 bytes and for- 
matted in ASCII code. 

Each record is headed with two lines of identification data. The first line contains 
the Marsden square number, the mean sub square (M SUB SQ) and the period of record 
(POR). Marsden squares are subdivided into 1 degree sub squares and this M SUB SQ shows 
the mean location of all of the observations within the analysis for this square. The POR is 
expressed as QORR-SSTT where QQ is the starting year, RR is the starting month, SS is the 
ending year and TT is the ending month. Generally, the POR is 7001-8412 indicating that 
the data is comprised of observations taken between January 1970 and December 1984. 

The second line contains the number of records examined for the analysis, the 
number of records accepted as valid, the number of records which failed the quality control 
checks and the number of records which failed the time constraints. 


Each subsequent line of data is preceded by a 10 character identification field and 
terminated with a6 character card number. The fonnat of the identification field is 
MSOMMDNT ET where MSQ is the Marsden square, MM is the month, DN is the day or 
night (1 for day. 2 for night) and TTT is the type of data. These data types are | for Paulus 
evaporation duct height, D1 2 for Pauius evaporation duct height crossed with Jeske duct 
height, W for wind speed. A for absolute humidity, M for modified air-sea temperature 
difference. r for rain rate. $ for gaseous attenuation rate at 35 GHz, t for gaseous attenua- 
tion rate at 94 GHz, u for total attenuation rate at 35 GHz, v for total attenuation rate at 
94 Gtiz, and M11 tor Paulus evaporation duct height crossed with wind speed. For the 
cross distributions, the MM does not represent the month but rather, represents an index 
field of the cross distribution. 

Details of the specific quantities and distribution characteristics are described in the 
following paragraphs. 


Paulus Evaporation Duct Height. The Paulus formulation of evaporation duct 
height calculations (reference 3) is a modification to the classical Jeske method (reference 4). 
The major difference in the approaches is that the Paulus technique attempts to account 
for inaccuracies in air temperature observations which are caused by thermal influence of 
the ship. 

The data are in month order subdivided by day and night. An observational set is 
drided inte 27 scetions (section number followed by data value). Sections |} through 22 ure 
the number of observations of evaporation duct height from 0 to 40 meters in 2 meter inter- 
vals with two additional categories: duct heights greater than 40 meters and duct heights 
that are not calculable (undefined). Sections 23 through 26 consists of the mean, first, 
second, and third quartile heights respectively. and section 27 is the total number of obser- 
valions during the period. 


Paulus Evaporation Duct Height Crossed with Jeske Duct Height. The joint 
probability of a duct height computed by the Paulus and Jeske.methods is represented by 
a 22 by 22 matrix where the 22 data line pairs and the first 22 intergers within each line 
pair represents height intervals from 0 to 40 meters in 2 meter intervals for both Paulus and 
Jeske duct heights. The real number and the final integer of each line pair represents the 
mean value of the quantity and the number of valid observations that the distribution is 
derived from respectively. 


Surface Wind Velocity. An observation set consists of 24 elements, i.e., the number 
of observations of wind velocity from O to 20 meters per second in | meter per second inter- 
vals; number of observations where winds exceeded 20 meters per second; the mean wind 
velocity for the period: and the total number of observations during the period. 


Absolute Humidity. An observation set consists of 24 elements, i.e., the number of 
observations of absolute humidity from less than 1 to more than 41 grams per cubic meter 
(g:m**3) in 2 grams per cubic meter increments: a mean absolute humidity for the period; 
and the total number of observations during the period. 


Modified Air-Sea Temperature Difference. The modified air-sea temperature differ- 
ence is defined in reference 3. Its purpose for inclusion within the DUCT63 database is to 
aid NOSC research efforts in the climatological description of evaporation duct heights. 

An observation set consists of 24 ele!ements, i-e., the number of observations of modified 
air-sea temperature difference from less than -10 to greater than +10 degrees Celsius in 

) degree increments; a mean modified air-sea temperature difference for the period: and the 
total number of observations during the period. 


Rain Rate. Rain rate is not a directly reported quantity in the STD-11 database. 
Rather, it is computed from the present weather code by techniques developed by Goroch 
(reference 5). An observation set consists of 24 elements, i.e., the number of observations 
of rain rate from less than J to greater than 5.01 mm/hr increments; the mean rate for the 
period; and the total number of observations during the period. 


Gaseous Attenuation Rate at 35 GHz. This quantity is derived from observations 
of air temperature, relative humidit: . and visibility. The reported visibility is used to 
calculate a value of the liquid water content based on the work of Johnson reported by 
Cook (reference 6). Methods described by Liebe (reference 7) are used to calculate the 
attenuation rate. Since the meteorological data is surface data, the attenuation rate is 
valid only for surface-to-surface propagation. It is not applicable for use with slant paths. 
An observation set consists of 24 elements, i.e., the number of observations of from less 
than 0.001 to greater than 5.01 dB’km in 0.025 dB/km increments; the mean gaseous 
attenuation rate for the period: and the total number of observations during the period. 


Gaseous Attenuation Rate at 94 GHz. The description above applies to the gaseous 
attenuation rate of 94 GHz with two exceptions. First the frequency is 94 GHz and second, 
the observations range from less than | dB/km to greater than 5.0] dB’km in 0.25 dB:km 
increments, 


Total Attenuation at 35 GHz. The major difference between the attenuation rate 
described for this quantity and the attenuation rate described in section 3 is that the total 
attenuation includes the effects of rain rate frorn section 3. The contributions from rain 
rate are modeled after the work of Falcone (reference 8). The data format is the same as for 
the 35 GHz gaseous attenuation rate. 


Total Attenuation at 94 GHz. The comments of the preceding section are directly 
applicable for the total attenuation rate at 94 GHz. Again, it should be stressed that tie 
meteorological data are observed at the surface and the attenuation rates described are 
strictly applicable to the surface-to-surface path. That is, the data are not applicable to a 
slant path. 


Paulus Duct Height Crossed with Wind Speed. The Paulus duct height crossed with 
wind speed is designed primarily for use by NOSC in the determination of evaporation 
ducting. The format of the data is identical to that as described within section 3 except the 
intergers within each line pair represents ‘vind velocities from O to 22 meters per second in 
1 meter per second intervals. 


SECTION 4 


HISTORICAL PROPAGATION CONDITION 
DATABASE CONSTRUCTION 


HEPC Radiosonde Observations 


By using the technique as described below, construction of the radiosonde based 
portion of the database is achieved. 


The following data (for an individual radiosonde station) are read from the GTE/ 
Sylvania Radiosonde Long A tape: 


Name Description data type 
GNSBD N-gradient first 100 meters nedian, monthly 
GNSTD N-gradient first 1OOOm_— rs M * 
PELOO “er occur elevated duct 002 % is 
PEL12 “e occur elevated duct 122 7 * 
PESOO “ occur elev-sfc duct 002 = “ 
PES12 “e occur elev-sfe duct 12Z ” * 
PSBOO “e occur sfe-base duct 00Z » He 
PSB12 “ occur sfc-base duct 122 Me 
AAG annual average gradient for elevated- average, annual 
surface and surface ducts in layer 

OHEL elev duct optm coupling hgt median, monthly 

MTEL elev duct median thickness se 7 

MDEL elev duct M-unit deficit ‘ 

GMEL elev duct M-unit gradient ” ie | 
MFEL elev duct trapping freq * a 
MTES elev-sfe duci thickness ” 7 | 
MDES elev-sfe duct M-unit deficit se = 
MFES clev-sfe duct trapping freq me i | 
MTSB sfe-buse duct thickness " ‘ 
MDSB sfe-base duct M-unit deficit i " 
MFSB sfc-base duct traping freq ~ + 

P2EL probability >1 elevated duct probability ” 

PSBEL probability sfe & elev duct . 1 

ACC00 number accepted soundings 00Z nuy oor” 

ACCI2 number accepted soundings 122 2 < 

NS median surface N-units median - 

SLEV station elevation (meters) 


AIL N-unit gradients are converted to M-unit gradients with the following convention: 
GMSBD = GNSBD + 156 (M ~ Km) (100 meter gradient) 

GMSTD = GNSTD + 156 (M- Km) (1000 meter gradient) 

TAG = (AAG +156) ; 1000 (M / Km) 


Figure 6 iNustrates the relationships of the above terms. 
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Figure 6. Variable definitions as used in elevated, surface-based and elevated- 
surface ducts. 


The elevated-surface and the surface-based duct data is combined to produce a 
surface-based duct with a positive lower gradient. This is accomplished by using a weighted 
average of the duct thickness, M-unit deficit and the trapping frequency as described below: 

#elev-sfe duct ESN = (PESOQ * ACCOO) + (PESI2 * ACC12) 

# sfc-based duct SBN = (PSBOO * ACCOO) + (PSBI2 * ACC12) 

# soundings with ducts TWD = ESN + SBN 


weighted thickness THK = (FSN * MTES + SBN * MTSB) / TWD 
weighted M deficit DM=-(ESN * MDES + SBN * MDSB) / TWD 
weighted trap freq T” ESN * MFES + SBN * MFSB) / TWD 


The optimum coupling heignt for surface ducts is computed as follows, insuring that 
a surface-based duct is formed from the statistics. 


optimum height OHSB = THK - DM/TAG 
if OHSB <0, then OHSB=THK / 2 


A surface-based duct M-unit deficit is computed as follows 
surface M deficit SMD =QOHSB * (GMSBD-1000) + DM 
If SMD is greater than zero, the weighted M deficit (DM) is recomputed such that 


the surtace deficit is 4 M-units. Setting the surface deficit to 4 M-units is a purely subjective 
decision. A surface-based duct must be created 


weighted M deficit DM =-(4 + SMD - DM) 


. Figure 7 illustrates the relationships of the weighted and adjusted statistical informa- 
tion tsed to construct a TESS acceptable surface-based duct. 


Table 2 illustrates the computed radiosonde database for World Meteorological 
Organization station 4YN, Fixed Ship, North Pacific Ocean. 
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YM, GMSBD (100 METERS) 


Figure 7. Variable definitions as used in 
M—o constructed surface-based duct. 


Table 2. Computed values for variables within the HEPC radiosonde database. 
WMO station 4YN, fixed ship North Pacific Ocean. 


Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Ann 


| 
SELV 12 12 12 12 12 12 12 12 12 12 12 12 12 
ACC daytime 116 115 129 #4117 #22132 «2118 75 107 106 95 108 106 1324 
ACC nighttime 129 112 127 «#107 «:133 ~=«120 &1 108 100 95 107 97 1307 
i NS 343 337 337 336 341 348 351 356 359 353 341 343 345 
4 GMSTD 107. 109 112 111 #109 4105 +=«105 98 102 104 108 107 107 
; GMSBD 85 89 92 71 97 97 94 72 &5 89 84 100 90 
OHSB 32 26 17 41 80 24 106 86 110 92 82 43 56 
THK 89 89 88 83 133 114 144 129 #137 #=121~=«#2125 58 98 
DM 4 4 5 7 11 6 14 10 13 12 10 8 9 
MFSB 730 960 746 577 575 1087 673 815 962 1014 692 3483 751 
PSB daytime 19 27 18 29 ll S 16 30 25 22 12 9 19 
PSB nighttime 5 3 7 9 2 1 1 3 1 3 3 0 3 
GMEL 128 122 #119 #119 #%116 #%&115 #4112 #%2112 2117 114 %121 #2120 «117 
OHEL 1286 1522 1561 1552 1522 1185 1515 1551 1572 1544 1542 1401 1498 
® MTEL 180 173 172 #4172 «+168 #203 210 220 171 #+4181 #215 #4177 «2186 
f MDEL | ea a ee ae a a ee a ee ek ee ee 
MFEL 206 186 201 205 197 151 133 119 212 216 #121 #4194 175 
PEL daytime 42 37 36 46 48 60 56 51 52 47 55 40 47 
PEL niazhttime 56 53 61 71 66 65 64 70 71 75 67 71 66 
P2EL 466 264 234 134 264 714 705 512 874 632 279 443 445 
PSBEL 339 264 195 848 226 294 641 791 348 368 186 99 376 
HEPC Surface Observation 


The HEPC function employs the Paulus evaporation duct heights and the surface 
wind velocities portion of the DUCT63. By using the technique as described below, 
construction of the surface observation database is achieved, 


: a. The following data (for an individual Marsden square} is read from the 
DUCT63 data tape: 


Name Description data type 


NHOD number height observations, day height, monthly 
NHON numbe. neight observations, night height, monthly 
TAOD total accepted observations, ducts monthly 


4 MSWV surface wind velocity ‘ean, mor niyv 
4 


NAOW total accepted observations, wind r-onthly 


16 


b. The percent occurrence of evaporation ducts is given by 


POED (day) =NHOD / TAOD 
POED (night) = NHON /TAOD 
ce. The mean surface wind velocity is transferred to the database without 
processing. 
d. If the total number of accepted observations is Jess than 10 for the evaporation 


duct height or the surface wind velocity, a -1. (indicating an error) is substituted for the 
calculated yalue. 


Tables 3 and 4 illustrate the computed surface observation database (day and night 
respectively) for Marsden square 85. 


Table 3. Computed values for variables within the HEPC surface observation database, 
Marsden square 85 (daytime). 


Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 
00-02 meters 3.1 2.4 2.2 1.7 2.3 2.8 3.2 3.6 1.9 2.2 1.8 1.8 
02-04 meters 1.8 1.8 1.8 8621.5 2.7 1.9 1.9 1.3 .7 1.2 21.6 
04-06 meters 5.2 4.6 2.9 2.7 2.7 4.9 5.6 3.9 3.6 2.0 3.1 3.4 
06-08 meters 7.6 7.6 6.3 4.7 5.9 7.6 7.6 7.2 5.7 5.4 4.3 7.0 
08-10 meters 11.0 12.2 11.0 9.6 10.2 11.4 12.4 10.8 9.8 8.9 8.1 10.3 
10-12 meters 14.5 15.6 14.9 12.8 17.6 16.8 19.2 15.3 13.3 12.7 12.6 13.7 
12-14 meters 16.6 16.1 17.4 15.8 18.3 19.0 18.4 17.6 17.0 16.0 15.4 15.0 
14-16 meters 13.7 13.4 13.8 15.1 15.8 13.6 14.7 16.6 16.4 16.5 16.1 14.5 
16-18 meters 11.5 10.4 12.7 14.8 11.6 11.1 8.5 11.4 12.9 12.2 12.5 12.1 
18-20 meters 7.2 6.9 7.8 10.6 7.4 5.5 4.4 5.6 7.810.3 10.9 8.9 
20-22 meters 4.4 4.5 5.2 6.5 4.1 2.5 2.4 3.6 4,7 7.1 7.0 5.6 
22-24 meters 1.9 2.5 2.3 2.9 1.7 1.2 1.1 1.6 3.1 3.0 3.8 3.4 
24-26 meters .8 1.0 1.3 9 6 4 3 6 1.4 1.4 1.3 1.7 
26-28 meters 4 5 4 5 .2 sd 1 aol" -5 1.0 1.0 .8 
28-30 meters .2 .2 2 3 a2 1 1 wl 13 4 .5 1 
30-32 meters .0 .2 20 .2 1 .0 0 .0 ot mae .2 sok 
32-34 meters .0 .0 .0 .0 0 .0 0 .O wl wl .0 wl 
34-36 meters 0 .0 0 .0 a) .0 0 .0 .0 .0 0 .0 
36-38 meters .0 0 0 .0 1 0 0 .0 .0 .0 .0 .0 
38-40 meters .0 0 0 .0 0 .0 0 .0 .0 me .O .0 
> 4O meters 6) 0 0 0 0 0 .0 .0 .O .0 .0 


mean height 1227 22.9: 133 ee! 1S ok P22 O12 6 19 oF NSS Ae 1d 7 
# of soundings 2278 2218 2397 2555 2480 2405 2355 2190 2074 2534 2391 2244 


wind (m/sec) 6.8 6.7 6.7 7.2 6.8 6.6 6.6 6.2 6.1 6.1 6.6 7.0 


Table 4. Computed values for variables within the HEPC surface observation database, 
Marsden square 85 (nighttime). 


Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 
00-02 meters 251. hed el 9° =#1.0 1.3 2.1 2.3 1.2 1.0 6 1.2 
02-04 meters 2.8 3.5 2.6 2.5 2.5 3.2 5.0 3.4 3.3 1.7 2.2 2.2 
04-06 meters 6.0 6.2 6.2 4.6 5.0 7.3 6.4 6.3 7.4 4.9 4.8 4.5 
06-08 meters 12.3 11.8 10.8 9.1 11.1 12.7 14.6 11.8 9.0 10.5 8.2 8.5 
08-10 meters 15.2 17.4 15.3 14.1 17.4 19.3 21.0 19.2 16.0 12.6 11.6 15.5 
10-12 meters 7.2 17.0 17.5 17.7 23.0 22.1 20.8 18.5 19.6 16.8 17.5 17.2 
12-14 meters 15.5 14.5 16.7 18.4 17.0 16.8 16.4 18.3 14.8 17.5 17.5 16.6 
14-16 meters 13.2 12.9 13.2 14.8 11.7 10°54 7.4 11.1 12.8 15.3 15.9 13.1 
16-18 meters 7.2 7.1 9.0 10.2 6.3 3.8 4.0 5.1 8.4 10.2 11.0 9.2 
18-20 meters 4.9 4.3 4.6 4.9 3.4 1.9 1.46 2.5 3.8 4.6 5.8 6.6 
20-22 meters 2.0 2.1 1.8 1.9 a) .9 6 1.0 2.4 2.8 2.4 3.3 
22-24 meters 8 ol 6 al i3 2 3 we 8 1.6 1.6 1.3 
24-26 meters 4 ae 5 ae i2 wl 2 .2 1 4 5 a) 
26-28 meters 2 wl 1 .2 .0 .0 0 lL 5 .3 a) 12 
28-30 meters 0 1 1 .0 .0 .0 0 .0 0 m4 ml wl 
30-32 meters 0 .0 0 .0 .Q .O 0 .0 0 wl 1 me 
32-34 meters 0 m0) 1 .0 .0 .0 0 Ll 0 oO 1 .0 
34-36 meters 0 .0 0 .0 .0 .0 0 .0 0 .0 .0 .0 
36-38 meters 0 .0 (6) .0 .0 .0 0 .0 0 .0 .O .0 
38-40 meters 0 .0 0 .0 .0 .0 0 .0 0 .O .0 .0 
> 40 meters 1 a) .0 0 0 0 .0 .0 0 .0 .0 wl 


mean height 11.4 11.3 11.7 12.0 11.3 10.5 10.1 10.7 11.5 12.3 12.6 12.3 
x of soundings 18/2 1750 1906 1970 2004 1927 1939 1789 1653 1998 1899 1770 


wind (m/sec) 6.5 6.5 6.5 6.8 6.5 6.4 6.3 5.9 5.7 5.7 6.2 6.7 


to 


La) 


6. 
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